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Abstract. We study general S 2 x S 1 Gowdy models with a regular past Cauchy 
horizon and prove that a second (future) Cauchy horizon exists, provided that a 
particular conserved quantity J is not zero. We derive an explicit expression for the 
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(3J^! metric form on the future Cauchy horizon in terms of the initial data on the past 

horizon and conclude the universal relation A p A{ = (8irJ) 2 where A p and A{ are the 
, areas of past and future Cauchy horizon respectively. 
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1. Introduction 



The well-known singularity theorems by Hawking and Penrose [20] show that 
cosmological solutions to the Einstein equations generally contain singularities. As 
discussed by Clarke [33] (see also [2] for a comprehensive overview) there are two types 
of singularities: (i) curvature singularities, for which components of the Riemann tensor 
or its kth derivatives are irregular (e.g. unbounded), and (ii) quasiregular singularities, 
which are associated with peculiarities in the topology of space-time (e.g. the vertex 
of a cone), although the local geometry is well behaved. In addition, the curvature 
singularities are divided up into scalar singularities (for which some curvature invariants 
are badly behaved) and nonscalar singularities (for which arbitrarily large or irregular 
tidal forces occur). The singularity theorems mentioned above provide, however, in 
general no information about the specific type of singularity — they make statements 
solely about causal geodesic incompleteness. This lack of knowledge concerning the 
specific nature of the singular structure is the reason for many open outstanding 
problems in general relativity, including the strong cosmic censorship conjecture and 
the BKL conjecture (see [1] for an overview). 

A major motivation for the study of Gowdy spacetimes as relatively simple, but 
non-trivial inhomogeneous cosmological models results from the desire to understand 
the mathematical and physical properties of such cosmological singularities. The Gowdy 
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cosmologies, first studied in [T7J [TS], are characterized by an Abelian isometry group 
£7(1) x £7(1) with spacelike group orbits, i.e. these spacetimes possess two associated 
spacelike and commuting Killing vector fields £ and 77. Moreover, the definition of Gowdy 
spacetimes includes that the twist constants e Q , J g 7 ,5£ a 7/ 3 V 7 £ l5 and e a ^s^ a r]^'V' r 'r] 5 (which 
are constant as a consequence of the field equations) are zerc0. 

For compact, connected, orientable and smooth three manifolds, the corresponding 
spatial topology must be either T 3 , S 3 , S 2 x S 1 or L(p, q), cf. |T8] (see also [26J 130| [15] ) . 
Note that the universal cover of the lens space L(p, q) is S 3 and hence this case needs 
not be treated separately, see references in [10] . 

In the T 3 -case, global existence in time with respect to the areal foliation time 
t was proved by Moncrief [25J. Moreover, he has shown that the trace of the second 
fundamental form blows up uniformly on the hypersurfaces t = constant in the limit 
t — > 0. As a consequence, the solutions do not permit a globally hyperbolic extension 
beyond the time t = 0. However, to date it has not been clarified whether the solutions 
are extendible (as non-globally hyperbolic C 2 -solutions) or are generically subject to 
curvature singularities at t — 0. 

Although global existence of solutions inside the "Gowdy square" (i.e. for < t < n, 
cf. Fig. [1] below) was shown by Chrusciel for S 2 x S 1 and S* 3 topology, see Thm. 6.3 
in [TO], it is still an open question whether globally hyperbolic extensions beyond the 
hypersurfaces t = or t = it exist. It is expected that these hypersurfaces contain either 
curvature singularities or Cauchy horizons; the theorem in [TO] however does not in fact 
exclude the possibility that these are merely coordinate singularities. 

For polarized Gowdy models, where the Killing vector fields can be chosen to be 
orthogonal everywhere, the nature of the singularities for all possible spatial topologies 
has been studied in [221 EI]- m particular, strong cosmic censorship and a version of 
the BKL conjecture have been proved. Investigations of singularities in the unpolarized 
case for T 3 topology can be found in [51 EH [321 [33] . 

For unpolarized S" 3 or S 2 x S 1 Gowdy spacetimes not many results on singularities 
(strong cosmic censorship, BKL conjecture, Gowdy spikes) are known. Particular 
singular solutions have been constructed with Fuchsian techniques in [34J. Moreover, 
numerical studies indicate that the behavior near singularities and the appearance of 
spikes are similar to the T 3 -case [161 El E] • 

In this paper, we study general (unpolarized or polarized) S 2 x S 1 Gowdy models 
with a regular Cauchy horizon (with S 2 x S 1 topology) at t — (cf. Fig. [TJ@ and assume 
that the spacetime is regular (precise regularity requirements are given below) at this 
horizon as well as in a neighborhood. As mentioned above, a theorem by Chrusciel 
[TO] implies then that the metric is regular for all t < ir, i.e. excluding only the future 
hypersurface t = it. With the methods utilized in this paper we are able to provide 

1 The assumption of vanishing twist constants is non-trivial only in the case of spatial T 3 topology. 
Note that in spatial S 3 or S 2 x S 1 topology there are specific axes on which one of the Killing vectors 
vanishes identically, which leads to vanishing twist constants. 

2 Without loss of generality we choose a past Cauchy horizon H p . 
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the missing piece, i.e. we prove that under our regularity assumptions the existence of 
a regular second (future) Cauchy horizon Hf (at t = n) is implied, provided that a 
particular conserved quantity J is not zero§. Moreover, we derive an explicit expression 
for the metric form on the future Cauchy horizon in terms of the initial data on the 
past horizon. From this explicit formula, the universal relation A P A{ = (8ttJ) 2 between 
the areas A p , A { of past and future Cauchy horizons and the above mentioned conserved 
quantity J can be concluded. 

The proofs of these statements can be found by relating any S 2 x S 1 Gowdy model 
to a corresponding axisymmetric and stationary black hole solution (with possibly non- 
pure vacuum exterior, e.g. with surrounding matter), considered between outer event 
and inner Cauchy horizon. Note that the region between these horizons is regular 
hyperbolic, i.e. the Einstein equations are hyperbolic PDEs in an appropriate gauge 
with coordinates adapted to the Killing vectors, see [21 [31 l23"fl (The Kerr metric is an 
explicitly known solution of these PDEs, see [31]. )@ As a consequence, the results on the 
regularity of the interior of such black holes and existence of regular Cauchy horizons 
inside the black holes obtained in [21 [31 [23] can be carried over to Gowdy spacetimes. 

The results in [2] were found by utilizing a particular soliton method — the so-called 
Backlund transformation. Making use of the theorem by Chrusciel mentioned earlier, it 
was possible to show that a regular Cauchy horizon inside the black hole always exists, 
provided that the angular momentum of the black hole does not vanish. (The above 
quantity J is the Gowdy counterpart of the angular momentum.) 

In [3l [23] these results have been generalized to the case in which an additional 
Maxwell field is considered. The corresponding technique, that is the inverse scattering 
method, again comes from soliton theory and permits the reconstruction of the field 
quantities along the entire boundary of the Gowdy square. Hereby, an associated linear 
matrix problem is analyzed, whose integrability conditions are equivalent to the non- 
linear field equations in axisymmetry and stationarity. Note that in this article we 
restrict ourselves to the pure Einstein case (without Maxwell field) and refer the reader 
to [31 [23] for results valid in full Einstein-Maxwell theory. 

We start by introducing appropriate coordinates, adapted to the description 
of regular axes and Cauchy horizons at the boundaries of the Gowdy square, see 
Sec. [21 Moreover, we revisit the complex Ernst formulation of the field equations and 
corresponding boundary conditions and introduce the conserved quantity J in question. 
In this formulation we can translate the results of [21 [31 E3] and obtain the metric on 

3 As we will see in Sec. I2.3[ the conserved quantity J vanishes in polarized Gowdy models. 

4 The interior of axisymmetric and stationary black hole solutions is non-compact and has spatial 
S 2 xE topology. Here the R-factor is generated by a subgroup of the symmetry group corresponding 
to one of the Killing fields. Therefore, it is possible to factor out a discrete subgroup such that S* 2 x S 1 
topology is achieved. 

5 Another interesting example of a spacetime with a region isometric to Kerr is the Chandrasekhar 
and Xanthopoulos solution 9 which describes colliding plane waves. It turns out that the region of 
interaction of the two waves is an alternative interpretation of a part of the Kerr spacetime region 
between event horizon and Cauchy horizon, cf. jT9j [21] . 
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the future Cauchy horizon in terms of initial data on the past horizon, see Sec. [3J As 
another consequence we arrive at the above equation relating A p ,A { and J, see Sec. HI 
Finally, in Sec. El we conclude with a discussion of our results. 



2. Coordinates and Einstein equations 

2.1. Coordinate system, Einstein equations and regularity requirements 

We introduce suitable coordinates and metric functions by adopting our notation from 
[j~6] . Accordingly, we write the Gowdy line element in the form 

ds 2 = e M (-dt 2 + d9 2 ) + sin t sin 9 [e L {d V + Qd5) 2 + e~ L d5 2 } , (1) 

where the metric functions M, L and Q depend on t and 9 alone. In these coordinates, 
the two Killing vectors are given by 

As mentioned in Sec. (TJ any S 2 x S 1 Gowdy model can be related to the spacetime 
portion between outer event and inner Cauchy horizon of an appropriate axisymmetric 
and stationary black hole solution. Black hole spacetimes of this kind have been studied 
by Carter [8] and Bardeen jl] . Among other issues they discussed conditions for regular 
horizons. In this paper we adopt their regularity arguments for our study of Gowdy 
spacetimes. Accordingly we rewrite the line element (P) in the form 

ds 2 = e M (-dt 2 + d9 2 ) + e 11 sin 2 fl(d^ + Qd5) 2 + e~ u sin 2 t d5 2 (3) 

where 

u = In sin t — In sin 9 + L. (4) 

Now, at a regular horizon (clear statements about the type of regularity follow below) the 
metric functions M, Q and u are regular, meaning that L possesses a specific irregular 
behavior therein 

At this point, some remarks about the specific regularity requirements needed in 
our investigation are necessary. A crucial role is played by a theorem of Chrusciel 
(Theorem 6.3 in [10]) which provides us with the essential regularity information valid 
in the interior of the Gowdy square. In this theorem it is assumed that initial data are 
given on an interior Cauchy slice, described by t = constant = to, < to < tt- These 
data are supposed to consist of (i) metric potentials that are if ^-functions of 6 and 
(ii) first time derivatives that are if fc_1 -functions of 9 (with k > 3). Here H k denotes 
the Sobolev space W k ' 2 that contains all functions for which both the function and its 
weak derivatives up to the order k are in L 2 . With these assumptions the theorem by 

6 We achieve the form of the line element used in [3J[3l[23] from ([3]) by introducing the Boyer-Lindquist- 
type coordinates (R,0,(p,t) with R :— recast, t := 5/(2^), = constant, and the metric functions 
fi := c M , u := c u , uj := — 2r\ L Q. Since the potentials jl > 0, u > and oj are regular at the axes and at 
the Cauchy horizon (cf. [4]), we see that M, u and Q are regular as well. 
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Figure 1. The Gowdy square. We assume a i? 4 -regular metric and £f 3 -regular time 
derivatives on all slices t = constant in a neighborhood (gray region) of the past Cauchy 
horizon (H p : t = 0) and hnd by virtue of the results in [2], that then the metric is ir- 
regular on all future slices t = constant, < t < 7r (unless the quantity J introduced 
in (1251) is zero). In particular, a _ff 4 -regular future Cauchy horizon (Tit : t = tt) exists. 



Chrusciel guarantees the existence of a unique continuation of the given initial data for 
which the metric is H k on all future spatial slices t = constant with to < t < tt, i.e. only 
the future boundary t = tt of the Gowdy square is excluded. (Note that Theorem 6.3 
as formulated in [10] assumes the metric to be smooth. However, this condition can be 
relaxed considerably to the assumption of H k spaces [12].) 

Now, for the applicability of our soliton methods it is essential that the metric 
potentials in ([3]) possess C 2 -regularity. Therefore, in order to apply both Chrusciel's 
theorem and the soliton methods, we need to require that the metric potentials M, u, 
Q be if 4 -functions and the time derivatives if 3 -functions of 9 on all slices t = constant 
in a neighborhood of the horizon H p , see Fig. [T]@ Then Chrusciel's theorem ensures the 
existence of an if 4 -regular continuation which implies (via Sobolev embeddings and the 
validity of the Einstein equations) that the metric potentials M, u, Q are C 2 -functions 
of t and 9 for (t, 9) G (0, tt) x [0, tt], i.e. in the entire Gowdy square with the exception 
of the two horizons H v {t = 0) and (t = n). Now, in accordance with Carter's 
and Bardeen's arguments concerning regularity at the horizon, we require that this C 2 - 
regularity holds also for t — 0, i.e. we assume in this manner a specifically regular past 
horizon % p . 

As mentioned above, these requirements allow us to utilize our soliton methods 
at Hp. Since 7f p is a degenerate boundary surface of the interior hyperbolic region, 
the study of the Einstein equations provides us with specific relations that permit the 
identification of an appropriate set of initial data of the hyperbolic problem at the past 
Cauchy horizon % v . 

7 In [21 [31 123) , the much stronger assumption was made that the metric functions be analytic in an 
exterior neighborhood of the black hole's event horizon. This stronger requirement was necessary to 
conclude that the metric is also regular (in fact analytic) in an interior vicinity of the event horizon, a 
requirement needed for applying Chrusciel's theorem. 
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For the line element ([3]), the Einstein equations read as follows: 

sin 2 

- u*t - cot t u t + u.ee + cot ua = 2 r^-e 2u (Q 2 t - Q 2 e ) , (5) 

sm t 

- Q,u + cot t Q tt + Q fie + 3 cot Q fi - 2(u )t Q )t - u fi Q fi ) = 0, (6) 

- M, tt + M,ee - ~u,t(u,t - 2 cot t) + \u fi {u, 9 + 2cot#) - \^^-e 2u (Q 2 - Q\) = 0.(7) 

2 2 2 sin i 

Alternatively to (J7J), the metric potential M can also be calculated from the first order 
field equations 

1 sin^ 

(cos 2 * - cos 2 0)M t = -e 2u [cos t sin 0(Q 2 + Q 2 ) - 2 sin t cos Q t Q e] 

2 sint L ' ' 

H — sin t sin [cos t sin ^(m 2 + u 2 e ) — 2 sin t cos # w^u^j 
+ (2 cos 2 1 cos 2 6* — cos 2 1 — cos 2 0)u j t 

+ 2 sin t cost sin ^ cos 9(u t e — tan. 9), (8) 
1 sin^ 9 

(cos 2 1 - cos 2 0)M 7 g = - -e 2u —j- [sin t cos 0(Q 2 + Q 2 e ) - 2 cos t sin Q, t Q, e ] 

— — sin t sin 6 1 [sin t cos ^(u 2 + u 2 e ) — 2 cos t sin UjUj] 
+ 2 sin t cos t sin cos # (m^ + tan t) 

+ (2 cos 2 1 cos 2 6» - cos 2 1 - cos 2 0)u e . (9) 



These expressions tell us that (see Appendix A for a detailed derivation) 

= Q,t — u,t — j Q — Qp — constant , M + u = constant (10) 

holds on H p . As the t-derivatives of all metric functions vanish identically at % p , a 
complete set of initial data at "H p consists of 

Q = Q p eR, ueH\ Q, tt eH 2 , (11) 

where Q t tt is in H 2 as a consequence of the regularity assumptions discussed above. 
Note that among the second t-derivatives only Q itt can be chosen freely since the values 
of M t ft as well as u# are then fixed, as again the study of the field equations ©-(Ej) 
near TL p reveals. Similarly, M is also fixed on H p by the choice of the data in (II ip . 

It turns out that the constant Q p is a gauge degree of freedom. This results from 
the fact that the line element ([T]) is invariant under the coordinate change 

E : (t, 0, ip, 5) £' : (t, 9,(p' = <p- tt5, 5), (12) 

leading to Q' p = Q p +£l in the new coordinates §. We use this freedom in order to exclude 
two specific values, namely Q p = and Q p = 1/J, where J is the already mentioned 
conserved quantity that will be introduced in fl26|) . This exclusion becomes necessary 
since the analysis carried out below breaks down if Q p takes one of these values. 

8 Note that for the corresponding black hole spacetimes, the coordinate change (|12[) describes a 
transformation into a rigidly rotating frame of reference (for more details see 
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We note further that as another consequence of our regularity requirements, the 
following axis condition holds at least in a neighborhood of the points A and B 
(cf. Fig. ED: 

A 1/2 : M — u. (13) 

Moreover, at these points A, B we have (see Appendix A) 

M A = M B = u A = u B . (14) 

Note that solutions which are also C 2 -regular up to and including % satisfy 
corresponding conditions at the points C and D. 



2.2. The Ernst equation 

In order to introduce the Ernst formulation of the Einstein equations, we define the 
complex Ernst potential 

£(jt,9)=f(jt,9)+ib(t,9), (15) 
where the real part / is given by 

/ := = -e~ u sin 2 t - Q 2 e u sin 2 # (16) 

and the imaginary part b is defined in terms of a potential a, 



a :-- 



CVi Q u ■ 2n 

— e sm V, 



via 



— sin t sin #60, 
f 



a« = — sin t sin 6 b t . 



;i7) 



(18) 



In this formulation, the vacuum Einstein equations are equivalent to the Ernst 
equation 

(-£, tt - cot t £ tt + E >ee + cot 9 E, ) = -£ 2 t + £ 2 e , (19) 

where $t(£) denotes the real part of £. As a consequence of (fl9l . the integrability 
condition a te = a^ t of the system (1TB1 is satisfied such that a may be calculated from 
f fTSj) using £. Moreover, given a and £ we can use f|T6|) and (fT7|) to obtain the metric 
functions w and Q. Finally, the potential M may be calculated from 

1 sin t sin 9 r . „ / „o -9 ,9 , 9 \ 
cos*sin0(/ 2 + /J + &5 + 6j) 



My 



/ 2/ 2 cos 2 t - cos 2 # 



-2sintcos^(/, t /, e + 6 t 6 e )-4/ 2 ^ 



cost 



sin t sin 



/> _ J 

/ 2f 2 cos 2 t - cos 2 9 

2costsm9(l t l e + b t b, e )-4f 2 



sin 9 

sint cos 9 (f 2 + f 2 +b 2 t + b 2 e ) 
cos9~ 



sint 



(20) 



(21) 



since the Ernst equation (|T9l also ensures the integrability condition M^e = M^ t . 
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As for the potentials introduced in Sec. I2.1l we conclude axis conditions which hold 
at least in a neighborhood of the points A and B (cf. Fig. [Q: 

Ay/2: £,o = 0, a = 0. (22) 

Moreover, at the points A, B we have / = 0. Again, solutions which are also if 4 -regular 
on "Hf satisfy corresponding conditions at the points C and D. 

It turns out that initial data £ p (9) = £(0,9) = f p (9) + ib p (9) of the Ernst potential 
are equivalent to the inital data set consisting of u, Q = Q p , Q^t at H p . Both sets are 
related via 

U = - Qy {m sm 2 ^, (23) 

r 9 

b p =b A + 2Q p (cos 9-l)-Q 2 l e 2u ^Q, tt (0, 9') sin 3 9' d9', (24) 

Jo 

where b A = b(0, 0) is an arbitrary integration constant. 



2.3. Conserved quantities 

As a consequence of the symmetries of the Gowdy metric, there exist conserved 
quantities, i.e. integrals with respect to 9 that are independent of the coordinate time 
t. One of them is J, defined by 

_ i r ojm e 2u( t , e) sin s e de = constant . (25 ) 

8 ,/ sint v ; 

As for the black hole angular momentum in the corresponding axisymmetric and 
stationary black hole spacetimes (cf. discussion at the end of Sec. H}, this quantity 
determines whether or not a regular future Cauchy horizon exists. In fact, it exists if 
and only if J ^ holds. Note that J vanishes in polarized Gowdy models, where we 
have Q tt = 0. 

It turns out that J can be read off directly from the Ernst potential and its second 
^-derivative at the points A and B on % p (see Fig. [I]), 

b B -^Q P ), Q P = --beeU (26) 
tt). 

A detailed derivation of these formulas can be found in |23|. 



where 

b B = b(t = 0,9 



3. Potentials on Ay, A2, and 

3.1. Ernst potential 

In the previous sections we have derived a formulation which permits the direct 
translation to the situation in which the hyperbolic region inside the event horizon 
of an axisymmetric and stationary black hole (with possibly non-pure vacuum exterior, 
e.g. with surrounding matter) is considered, as was done in [2"l [31 12"3]. 
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In [2] it has been demonstrated that a specific soliton method (the Backlund 
transformation, see Appendix B) can be used to write the Ernst potential S in terms of 
another Ernst potential So which corresponds to a spacetime without a black hole, 
but with a completely regular central vacuum region. Interestingly, the potential 
So = So(t,6) possesses specific symmetry conditions which translate here into 

£o(t,Q) = £o(0, t) potential at A\, 

So(t, 7r) = £o(0, it — t) potential at A2, 

So(n, 6) = £o(0, n — 9) potential at Hf. 

Hence the potential values at the boundaries Ai, A2 and Hf are given explicitly in terms 
of those at H p . Now the Backlund transformation carries these dependencies over to the 
corresponding original Ernst potential S, i.e. we obtain S at A±, A2 and Hf completely 
in terms of the initial data at H p . 

An alternative approach (see [3j [23]) uses the inverse scattering method. In these 
papers the potentials on A\, A2 and Hf were obtained from the investigation of an 
associated linear matrix problem. The integrability conditions of this matrix problem 
are equivalent to the non-linear field equations, see Appendix B We may carry the 
corresponding procedure over to our considerations of Gowdy spacetimes. Accordingly 
we are able to perform an explicit integration of the linear problem along the boundaries 
of the Gowdy square. Since the resulting solution is closely related to the Ernst potential, 
it provides us with the desired expressions between the metric quantities on the four 
boundaries of the Gowdy square. 

Note that in both approaches the axes Ai and A2 are considered first. Starting at 
H p and using the theorem by Chrusciel |10j . which ensures if 4 - regularity of the metric 
inside the Gowdy square (i.e. excluding only Hf), we derive first the Ernst potentials 
at A\ and A2 in terms of the values at H p . It turns out that for J 7^ these formulas 
can be extended continuously to the points C and D at which A\ and A2 meet Hf 
(cf. Fig. [1]). Moreover, with the values at C and D it is possible to proceed to Hf, and 
in this way we eventually find an Ernst potential which is continuous along the entire 
boundary of the Gowdy square. As the theorem by Chrusciel ensures unique solvability 
of the Einstein equations inside the Gowdy square, we conclude that the if 4 -regularity 
of the Ernst potential holds up to and including Hf which therefore turns out to be an 
if 4 -regular future Cauchy horizon. 

The resulting expressions of the Ernst potentials at the boundaries A\, A2 and Hf 

read 

A p ( \ eu a n\ i[&4 - 2Q P (x - l)]S p (x) + b\ 

*(* = «ccosM = 0) = S p (x)-i[b A + 2Q p ( X -l)] ' (27) 

A 2 : S 2 ( X ) := Sit = anxttsC-s), e = *) = % B r" Q ttl £ fl + ^ » ^ 

S p (x) -i[b B + 2Q p {x + 1)J 

~/ s ~, „ , ai{x)SJx) + a 2 (x) , , 

Hf : Sf{x) := S{t = 71,6 = arccos -x = 1 p 21 ' 29 

bi( X )S v {x) + b 2 {x) 
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where 

S p (x) := £(t = 0, 9 = arccosx) (30) 

denotes the Ernst potential on "H p and a±, 02, 61, and b 2 in (l29i) are polynomials in x, 
defined by 

ai = i[16Q 2 (l - x 2 ) + 8Q p {b A {x + 1) + b B {x - 1)) 

+ (6a " b B )(b A (x - l) 2 - + l) 2 )], (31) 

a 2 = 8Q p [b 2 A (x + 1) + 6|(a; - 1)] - Ab A b B (b A - b B )x, (32) 

61 = 4(4Q P + 6 B - 6 A )x, (33) 

b 2 = i[4Q p (l - x 2 ) - 6 A (1 + x) 2 + 6 B (1 - a;) 2 ](4Q p + b B - b A ). (34) 

A discussion of (1291) shows that £f is indeed always regular provided that the black 
hole angular momentum does not vanish, which in turn means that J 7^ 0, cf. (j25i) . In 
order to prove this statement, we first note that both numerator and denominator on 
the right hand side of (1291) are completely regular functions in terms of x, since ai, a 2 , &i, 
62 are polynomials in x and the initial function S p is regular by assumption. Hence, an 
irregular behavior of the potential £f could only be caused by a zero of the denominator. 
Consequently, we investigate whether the equation 

h(x)£ p (x) + b 2 (x) = (35) 

has solutions x 6 [—1, 1]. The real part of fl35l) is given by 

4a;(4Q p + b B - 6 A )/ P (x) = 0. (36) 

Using ( 1231) and ( 1261) together with our gauge Q p 7^ and the assumption J^Owe find 
that (|36l) has exactly the three zeros, x = —1, x = and x = 1 (corresponding to = 71, 
9 = 7r/2 and = 0). Now, for a; = the imaginary part of ( 1351) does not vanish, whereas 
for x = ±1 it does. Thus we find that the only zeros of the denominator in (129]) are 
located at the two axes (x = ±1). As a matter of fact, the regular numerator of (1291) 
also vanishes at x — ±1, as can be derived in a similar manner. Consequently, we study 
the behavior of Sf for x = ±1 in terms of the rule by L'Hopital. As both numerator 
and denominator in (1291) have non-vanishing values of the derivative with respect to x 
for x = ±1, we conclude that the Ernst potential is regular everywhere whenever J 7^ 
holds. 

Consider now the limit J — > for which the expression 
4Q P + b B -b A 

vanishes, cf. (1261) . As this term appears as a factor in both b\ and b 2 (cf. (|33l) . (l34l) ). we 
find that the denominator in (|29|) vanishes identically. The numerator, however, remains 
non-zero, which means that Ernst potential diverges on the entire future boundary t = 7r, 
< 9 < 71. We conclude that "Hf becomes singular in the limit J — > 0. This divergent 
behavior of the Ernst potential corresponds to the formation of a (scalar) curvature 
singularity at H{. In order to illustrate this property, we calculate the Kretschmann 
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scalar at the point C on "H f (see Fig. [T]). Using the axis conditions discussed in Sec. [2] 
and the Einstein equations, we obtain 

Ri jk iR m \c = 12 [e" 2u (l + 2u tt ) 2 - Q%] c . (37) 

In terms of the Ernst potential, this expression reads (cf. Eq. (1401) below) 

Ri jk iR i3kl \c = I [(f,tt + f,uuf - (btt + butt) 2 } c ■ (38) 

Now we can use (1271) to derive a formula that contains only the initial data on the past 
horizon 7-L p . Together with ( 1261) we get 

Ri jM R ijkl \c = -^T/e [^QP 2 - ^Ql J ~ & - ^QpV,** - V,*Y] B , (39) 

where x = cos 8. Note that the numerator is well-defined and bounded for our ir- 
regular metric, a fact which is ensured by the validity of the Einstein equations near 
T-i p . 

Equation (|39|) indicates that the Kretschmann scalar diverges as J -6 in the limit 
J — > 0. In fact, as we choose Q p ^ (see Sec. I2.ip . and furthermore f tX ^ holds 
(because 2nf jX = —Q 2 A P where A p , < A p < oo, is the horizon area of 7i p , see Sec. Hj), 
we conclude that f\ is the dominating term in the numerator of (139]) for sufficiently 
small J. Hence the Kretschmann scalar indeed diverges as J~ 6 in the limit J — > 0. 



Metric potentials 

From the Ernst potentials £ x = h + i6 x , £ 2 = f 2 + ib 2 , £r = ft + ib f in (f2Tj) , (l28]h (|29|) 
we may calculate the metric potentials M, Q and u on the boundaries of the Gowdy 
square. Using (HI, (HZ]), (HD, (HI, (PD, (H we obtain 

^ : e Ml = e- = Q 1 = (40) 



A:e^=e- = -^, Q 2 = (41) 
/2 2 sin t 

% = e* = -^^, Q = Q f , e- = -— * (42) 



where 



b A -b B - 4Q P 
Note that Qf 7^ in our gauge (cf. (T2"6"D ): 

&a - &b + 4Q P = (6 A - b B - 4Q P ) + 8Q P = -%Q p J + Q P = 8Q P (1 - JQ P ) ^ 
in accordance with the discussion in Sec. 12. II where the gauge freedom was used to assure 
7^ Q P 7^ I /J. Furthermore, using ( I2"7j) - (I291) and our regularity assumptions for the 
initial data, it is straightforward to show that (— sin 2 t//i), (— sm 2 t/ f 2 ) and (— sin 2 9/ ff) 
are regular and positive functions on the entire boundaries Ai, A 2 and Hf, respectively. 
Moreover, the terms (&i,t/ smt) and {b 2it / sint) are regular on A\ and A 2 , respectively. 
Consequently, the above boundary values for the metric potentials M, u and Q are 
regular, too. 
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4. A universal formula for the horizon areas 



In [2] a relation between the black hole angular momentum and the two horizon areas 
of the outer event and inner Cauchy horizons was found. This relation emerged from 
the explicit expressions of the inner Cauchy horizon potentials in terms of those at the 
event horizon. Translated to the case of general S 2 x S 1 Gowdy spacetimes, this relation 
is given by 

A P A{ = (8vr J) 2 , (44) 

where the areas A p and Af of the Cauchy horizons % v and Hf are defined as integrals 
over the horizons (in a slice 5 = constant), 

TT 

A p /{ = J y/goo9w d#d^ = 2vr J e^ | Wp/f sin 6d6 = 4ire u \ A /c- (45) 
s 2 o 



5. Discussion 



In this paper we have analyzed general S 2 x S 1 Gowdy models with a past Cauchy 
horizon % p . As any such spacetime can be related to a corresponding axisymmetric 
and stationary black hole solution, considered between outer event and inner Cauchy 
horizons, the results on the regularity of the interior of such black holes (obtained in 
[21 El [23]) can be carried over to the Gowdy spacetimes treated here. In particular, spe- 
cific soliton methods have proved to be useful, (i) the Backlund transformation and (ii) 
the inverse scattering method. Both methods imply explicit expressions for the metric 
potentials on the boundaries A\, A2, Hi of the Gowdy square in terms of the initial 
values at H p . Moreover we obtain statements on existence and regularity of a future 
Cauchy horizon as well as a universal relation for the horizon areas. These results are 
summarized in the following. 

Theorem 1. Consider an S 2 x S 1 Gowdy spacetime with a past Cauchy horizon 7i p , 
where the metric potentials M, u and Q appearing in the line element (Ej) are H A - 
functions and the time derivatives H 3 -functions of the adapted coordinate 9 on all slices 
t = constant in a closed neighborhood N := [0, to] x [0, tt], to 6 (0, tt), of Hp. In addition, 
suppose M,Q,u £ C 2 (N). Then this spacetime possesses an H 4 -regular future Cauchy 
horizon T-L{ if and only if the conserved quantity J (cf. < 0?3)) ) does not vanish. In the limit 
J — > 0, the future Cauchy horizon transforms into a curvature singularity. Moreover, 
for J 7^ the universal relation 

A P A { = (8?r J) 2 (46) 

holds, where A p and Af denote the areas of past and future Cauchy horizons. 

Remark. Note that the statements in Thm. 1 can be generalized to 5* 2 x S* 1 Gowdy 
spacetimes with additional electromagnetic fields, see [31 [23] . The proof utilizes a more 
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general linear matrix problem in which the Maxwell field is incorporated. Again the 
corresponding integrability conditions are equivalent to the coupled system of field equa- 
tions that describe the Einstein-Maxwell field in electrovacuum with two Killing vectors 
(associated to the two Gowdy symmetries). It turns out that apart from J a second 
conserved quantity Q becomes relevant. The corresponding counterpart of this quan- 
tity in Einstein-Maxwell black hole spacetimes describes the electric charge of the black 
hole. For Gowdy spacetimes we conclude that a regular future Cauchy horizon exists 
if and only if J and Q do not vanish simultaneously. Moreover, we find that Eq. ( 1461 ) 
generalizes to A P A { = (8ttJ) 2 + (AttQ 2 ) 2 . 

With the above theorem we provide a long outstanding result on the existence of a 
regular future Cauchy horizon in S 2 x S 1 Gowdy spacetimes. We note that the soliton 
methods being utilized in order to derive our conclusions are not widely used in previous 
studies of this kind. Therefore we believe that these techniques might enhance further 
investigations in the realm of Gowdy cosmologies. 
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Appendix A. Derivation of initial and boundary conditions 

We provide a derivation of the initial and boundary conditions through a thorough study 
of the Einstein equations near the past Cauchy horizon % p , i.e. the initial surface t = 0. 
First multiply the field equation ([6]) with sint and consider subsequently the limit t — > 0. 
Taking our regularity assumptions into account (cf. discusscion in Sec. I2.ip . we arrive 
at 

Q t = at t = 0. (A.l) 

With the result ( 1A.1I) . we study the limit £ — >• of (jEj) in terms of the rule by L'Hopital 
and find 

sin 3 9 e 2u Q i9 = constant. (A. 2) 

Evaluation at 9 = shows that the constant vanishes, leading to 

Q = constant at t = 0. (A. 3) 

Next multiply Eq. (|5]) with sint and study the limit t — > 0. With (1 A. II) and (1A.3I) we 
obtain 



u t = for t = 0. 



(A.4) 
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(A.5) 



Using the previous results, we derive from OH]) in the limit t — > 

Af t = at t = 0. 
On the other hand, @ leads to 

M + u = constant at t = 0. (A. 6) 

Similarly, we study the Einstein equations on the axis, i.e. in the limit sin6 — > 0. 
Multiplication of Eqs. (E]) and ([6]) with sin 9 leads for sin 9 — > to 

w e = 0, Q je = for sin# = 0. (A.7) 

From OH]) and (Q we obtain 

Mo = 0, M-u = constant for sin0 = O. (A.8) 

As a consequence of an axis regularity condition that excludes the appearance of struts 
or knots along the axis (see [35] for details), it turns out that the constant in (1 A. 81) 
vanishes. Hence we get 

M = u for sin# = 0. (A.9) 

At the points A and B (see Fig. [1]), (1A.9j) means that Ma = ua and Mb = ub, and 
with Ma + ua = Mb + ub (cf. ( 1A.6I) ) we conclude 

Ma — ua — Mb — Ub- (A.10) 
Finally, we derive 

6,e = for sinfl = (A.ll) 

by multiplying the Ernst equation f|T9]) with sin^ and considering the limit sin# — > 0. 
Moreover, it follows from the definition ( ITTj) of the potential a that 

a = for sinfl = 0. (A. 12) 



Appendix B. The linear problem and Backlund transformations 



In this appendix we briefly discuss the mathematical structure of the Ernst equation (TT9|) 
which permits the application of so-called soliton methods. More details can be found 
in [2] |3] [23]. For a sophisticated introduction to soliton methods for the axisymmetric 
and stationary Einstein equations we refer the reader to [29J. 

There are two soliton methods which lie at the heart of the treatment of S 2 x S 1 
Gowdy spacetimes pursued in this paper: (i) the Backlund transformation and (ii) the 
inverse scattering method. Both methods make use of the following linear matrix 
problem (see [27] [28]), which read in our coordinates as follows: 



B x 





B,„ 



A x 




+ A 



B x 
A~ 



B„ 



A A 







(B.1) 
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Here, <fr = <&(x,y,K) is a 2 x 2 matrix pseudopotential depending on the coordinates 
x = cos(t-0), y = cos(t + 9) (B.2) 
as well as on the spectral parameter K G C. The function A is defined as 



A(a;,y,JO = y§rf- (B.3) 

For fixed values x, y, the equation (IB. 31) describes a mapping C — > C, >-)■ A from a 
two-sheeted Riemann surface (i^-plane) onto the complex A-plane. In the i^-plane the 
two .fT-sheets are connected at the branch points 

K 1=x (A = oo), K 2 = y (A = 0). (B.4) 

Examining the integrability conditions <&, X y — *&,yx yields, on the one hand, that the 
quantities A x , A y , B x and B y are given in terms of a single complex 'Ernst' potential 
£ = f + ib, 

On the other hand, the integrability conditions xy = *&,yx tell us that this potential 
£ satisfies the Ernst equation ( lT9i) . Conversely, any solution £ to the Ernst equation 
implies the existence of an associated matrix $ which obeys the above linear matrix 
equations (IB. If) where the functions A x , A y , B x and B y follow from (IB. 51) . 

Now, with a Backlund transformation a new potential £ can be constructed from a 
previously known one Sq. Starting from £ and the corresponding matrix function <fr , 
we consider transformations of the form 

BT n : $ ^ * = T n $ , n e IN even, (B.6) 

where T n is a matrix polynomial in A of degree n. From determined via (IB. 61) . one 
can finally calculate the corresponding new Ernst potential £, see [29] . 

Note that our specific expressions for the metric at the future Cauchy horizon Hf 
in Gowdy spacetimes can be obtained by considering the particular case of a twofold 
Backlund transformation (n = 2), for which the new Ernst potential £ reads 

g [ati (cos t + cos 9) — a 2 (cos t — cos 9)} £ + 2£ 
aii ( cos t + cos 9) — a 2 (cos t — cos 9) — 2 
Here, ot\ and «2 are solutions of the Riccati equations 

<V= - (Aiaf + + (ai + A i )-^ ) (B.8) 



1 2 \ &0,y ( 1 \ £o,y 



«i*= -(-«?+«*) + r i + AJ 2/b' ' = 1,2 ' (R9) 



with 
where 



Ai := X(x,y,K = -1), A 2 := X(x,y,K = 1). 
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In our second approach, the inverse scattering method, the linear problem (IB. 1[) 
is integrated along the boundaries of the Gowdy square. It turns out that explicit 
formulas can be found and that, moreover, the resulting solution must be continuous at 
this boundary (provided that the solution is regular at "Hf, which is true for J ^ 0, see 
discussion in Sec. [3]). In this way we find the expressions that constitute the statements 
of this paper. 
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